Improving quantum entanglement through single-qubit operations 
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We show that the entanglement of a 2 x 2 bipartite state can be improved and maximized prob- 
abilistically through single-qubit operations only. An experiment is proposed and it is numerically 
simulated. 



PACS numbers: 03.65.Ud, 03.67.Ac 

Introduction. — Quantum entanglement plays a central 
role in quantum information and also in the foundations 
of quantum physics. Thus, it has been extensively stud- 
ied (see, e.g., [l|-|6|). One important topic here is how to 
improve quantum entanglement of a bipartite quantum 
state As is well known, quantum entanglement can 
be improved through entanglement purification Q where 
a bipartite state is first transformed to a Werner state 
and then two-qubit operations at each sides are needed 
to improve the quantum entanglement probabilistically. 

In this letter, we shall present a theorem (Theorem 2) 
to maximize the entanglement of a two-qubit mixed state 
through single-qubit operations only. The theorem can 
be used to efficiently improve the quantum entanglement 
of a mixed state without the difficult 2-qubit operations. 
Explicitly, given a two-qubit mixed state p\ n = pyi, by 
taking local (non-trace-preserving Q) maps on qubit 1 
and qubit 2 separately, what is the maximally achievable 
entanglement of the normalized outcome state, and what 
are the specific maps needed on each qubits. 

To make a clear picture of our work we consider the 
following example with a pure state p- ln = \x){x\ and 
|x) = a|00) + b\li) first. Take the following specific non- 
trace- preserving map on the first qubit: 

e ® I(p- m ) = M{h, b) <g> I ■ p m ■ M (5, b) <g> I, (1) 

■6|1)(1|, and |a| 2 + H 2 = 1. We 



where M(a, b) = a 
have 



e®I(pin) = l\x')(x'\ 



(2) 



and 7 



\bb\ 2 , \x') 



|00) + ^|11). The 



entanglement concurrence of the outcome state is 



c(\x')ix'\) 



2\aabb\ 



\bbf 



(3) 



Setting \a\ = |6| and \b\ = \a\, we shall obtain the max- 
imum output entangled state \cj) + ) = -^(|00) + |11)) 
state (up to a normalization factor). Physically, the map 
M(a,b) can be easily realized. For example [Pj, one can 



use a polarization-dependent attenuator, with transmit- 
tance proportional to a for a horizontally polarized pho- 
ton (state |0)) and transmittance proportional to b for a 
vertically polarized photon (state |1)). Once we find a 
photon at the outcome port of the attenuator, the ini- 
tial state |x)(xl has been mapped to the outcome state 

Ix'Xx'l- 

Outline of our work. — Our goal is to look for the largest 
achievable entanglement through local operations, i.e., 
among all physical maps e® e' : p ou t — £ ® 6 '{pin)j which 
map gives out the largest entanglement of the outcome 
state pout ■ Most generally, any local map e <g> e' can be 
represented in the form of Kruss operators Q : 



Pout = r * ® r '> • • ( r * r '>) f = J2 p *P i 



(4) 



where ptpi = Ti (8 r-pi n (Fj <g) T-) 1 '. Denote C{p) as the 
entanglement concurrence of state p. Suppose C(p m ) 
is the largest among all {C(pi)}. Obviously, 



C(pout) < $>iC(pi) < C( Pm ). 



(5) 



Therefore, to find the largest entanglement concurrence 
of the outcome state among all local maps, we only 
need to seek it in the following special class of maps: 
Pout = Q <E> Q'pin(Q ® Q') f wh ere Q, Q! are 2 x 2 posi- 
tive matrices. According to the singular-value decompo- 
sition, the positive matrix Q(Q') can be decomposed into 
Q = UxDUiQ' = U[D'U') where U^U^U') are uni- 
tary matrices and D(D') is a positive-definite diagonal 
matrix. Since a unitary transformation plays no role in 
the entanglement, we only need to consider the positive 
matrices in the form of DU (D'U 1 ). 

As shown in Lemma 2, any two-qubit state p ln can be 
generated from the maximally-entangled state \cf> + ) acted 
by a one-sided map I <g> e' , i.e., p m = I (8 e'{\<j) + ) {cj) + \). 
Therefore, we can start with entanglement evolution and 
maximization under non-trace-preserving one-sided maps 
and then apply the result to the general problem of im- 
proving and maximizing quantum entanglement through 
single-qubit operations. 
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Entanglement evolution and maximization under non- 
trace-preserving maps. — A non-trace-prcscrving one- 
sided map I <g> d is fully characterized by p e < = I ® 
e'd^+X^+Q [Tl| . We assume 

I®e'(\<t>+)(<f>+\) = fp E ,, /®e'(|V)(Vl)=/W . (6) 

where / = tr [I <g> £'(|^ + )(^+|)], /' = tr [I <g> e / (|^)(^|)]. 

A 2 x 2 pure state |x) = a|00) + &[ IX) can be rewritten 
in the form |x> <xl = 2M(a, b) <g> I(\4> + ){(j>+\)M{a, b) ® 7. 
From Eq. © we have 7 <g> e'(|x)(x|) = 2fM(a,b) ® 
Ip e 'M(a,b) ® I = f'p x - We emphasize here that 
even though p e / is normalized, the operator 2M (a, 6) ® 
Ip e /M(a, b) ® 1 is not necessarily normalized. Define the 
following function C of an arbitrary non- negative definite 
4x4 matrix (operator) N 

C(N) = max{0, A ~ V 7 ^ ~ V^ 3 ~ V^J, (?) 

where are the eigenvalues of JV -N, in descending or- 
der, with N = a y ® a y N*a y ® a y , and N* is the complex 
conjugate of N. If TV is a density matrix of a 2 x 2 sys- 
tem, C(N) is just the entanglement concurrence of the 
system [10(. With this definition of C, we can summarize 
the major result, equation (5) in Ref. [l2[ as: 
Lemma 1. Given any density matrix p e >, if N = 
2M(a, b) ® Ipe'M^a, b) ® I, then 



C(N) = C(\ X )( X \)-C(p £ 



2\ab\C(p £ 



(8) 



However, this is not the entanglement concurrence of 
p x because N is not necessarily normalized, even though 
p s i is. Now denote N = gp x , and g = triV. According to 
the definition of C and p^ in Eq. (JB]), 



C{p x ) = C(N)/g = 2\ab\C(p £ ,)/g 



(9) 



where g = tiN = 2tr[M(a, b)®Ip £ ,M^(a, b)®I]. To avoid 
meaningless results, we assume C{p' £ ) > throughout 
this paper. Assume that the density matrix of the first 

°\ a I = Kn, where tr 2 
a* 01 ) 

is the partial trace over the subspace of the second qubit 
and Ci = (0|tr2/v|0), c 2 — (l|tr2p e '|l). Consequently, 



qubit of p £ i is po = tr 2 p e 



g = 2tr[M (a, 6)/? M + (a, 6)] = 2|a| 2 Cl + 2|&| 2 c 2 . (10) 
Therefore, the value of output entanglement 

2\ab\C{p £ ,) 



C[Px) - (|ap Cl + |&Pc 2 )' 



is maximized when \a\ = 
mum value 



Ea, |6| 



(11) 

cT, with the maxi- 
(12) 



2 y ?C 1 C 2 

More generally, the initial pure state can be 

|V) =1® U\ X ) = V2M(a,b)®U\4>), (13) 



where U is an arbitrary unitary operator. Given the fact 
that U* ® U\<fi + ) = \(f> + ) for any unitary U, we have 

V2M(a,b)®U\(j) + ) = \/2M(a,b)U T ®I\4>+). (14) 



In such a case, we obtain 



C{p^)=2\ab\-C{p £ ,)/g' 



(15) 



and g' = tv[M{a,b)U T ® Ip s >U*M(a,b) ® I}. To max- 
imize C(p^), we first fix U and maximize it with a, b. 



Assume U t KqU* 



The largest value for 



C(p^) is C(p £ ')/2 x /c' 1 c' 2 , as shown already. To maximize 
the value over all U, we only need to minimize d 1 d 2 - Since 
U is unitary, det(U T K U*) = det K . Therefore d t d 2 = 
det Kq + |c/| 2 , which is minimized when a' = 0. Namely, 
C(p^) is maximized when U t KqU* is diagonalized and 
y/k = c' 2 , Vb = ci, i.e., M(a,b)U T (tr 2 /v) U*M{a,b) = 
diag[l/2, 1/2]. We obtain 

Theorem 1. Denote Q to be a 2 x 2 positive-definite 
matrix. Given the inseparable two-qubit density matrix 
Pin = Pe' = I ® £ ' (\(f> + )(4> + \), the entanglement of the 
normalized density matrix p\ = Q® Ip m {Q ® Iy maxi- 
mizes when Q (tr 2 p in ) = diag[l/2, 1/2] and the entan- 
glement concurrence is: 



C 



M 



C(p £ >) 



2^det [tr 2 p £ 



(16) 



Improving and maximizing quantum entanglement 
through single-qubit operations. — To apply our theorem, 
we need the following lemma: 

Lemma 2. Given any 2x2 bipartite mixed state P12, 
there exists a map d such that p ln = I ® d ' (\<fi + }((j) + \). 

Note that map d here is in general non-trace- 
preserving. Since any two-qubit density matrix p- ln can 
be decomposed into the mixture of a few pure states, 
say pi n = Aj|^i)(V'i|- Obviously, for any bipartite 
pure state there always exists a positive opera- 

tor M[ such that |^) = I ® M-|0 + ). Therefore, we 
have p in = Ei^i 7 ® Ml\(f) + )(<j) + \I ® M,? f . Denoting 
/®e'(|^+>(0+|) = Ei 7 ® VAiM/l^^+l/® VAiMf 
completes the proof. 

With Theorem 1 and Lemma 2, we can improve the 
quantum entanglement of any 2-qubit state p- ln (here and 
after, the 2-qibit states are normalized) step by step, 
with single-qubit operations only. Denote K\ = tr2Pi n , 
if det K\ < 1/4, we construct Mi{a\,b\) and local uni- 
tary Ui such that M^KxUImI = diag[l/2, 1/2]. The 
local operation on qubit 1 transforms state p- m into 
the outcome state p\ = M\Ui ® Ip lr Jj\Ml ® I. Ac- 
cording to Theorem 1, the entanglement concurrence 
of the outcome state is C{p\) = C(pi n )/(2Vdet K{) > 
C(pi n ). The normalized density matrix of qubit 1 is 
tr 2 pi = diag[l/2, 1/2] now, but in general the normal- 
ized density operator of qubit 2 is not diag[l/2, 1/2] 
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now. Using Lemma 2, we know p\ can be written 
in the form of pi = e ® I(\4> + )(cf> + \). We can now 
apply Theorem 1 again to further improve the quan- 
tum entanglement through operation on qubit 2. De- 
note K[ = tripi. If det(K[) < 1/4, we construct 
new operators M[ and U[ such that the density ma- 
trix of qubit 2 is diag[l/2, 1/2] after the operation, i.e., 
M[U[{tT: 1 p 1 )u' l Hl^ = diag[l/2,l/2]. The operation on 
qubit 2 leads to a new outcome state p[ = I<£)M[U[piI® 
Ui . The operation on qubit 2 improves the entan- 
glement concurrence to C(p' 1 ) = C(pi)/(2^/dct K^) > 
C(pi). After the non-trace-preserving operation above 
on qubit 2, we have K[ = tiip[ — diag[l/2, 1/2], but in 
general the density matrix of qubit 1 is not diag[l/2, 1/2], 
i.e. K2 = t?2p'i 7^ diag[l/2, 1/2] now. We can construct 
new operators Mi and U2 to improve the entanglement 
of p' x . The process will continue step by step until the 
determinant of two reduced density matrices are all equal 
to 1/4 after many steps of iterations. Since the entan- 
glement concurrence of a two-qubit state can never be 
greater than 1 and the entanglement always increases 
during the iteration process above, there must exist a 
limit value of the entanglement in the process say, after 
many steps of iterations, the process gives out the largest 
entanglement concurrence. This also means that after 
many steps of iterations, the process always produces a 
two-qubit state where the reduced density matrices of 
each qubit are diag[l/2, 1/2] simultaneously. Therefore 
the process that transfers the reduced density matrix of 
qubit 1 and the reduced density matrix of qubit 2 into 
diag[l/2, 1/2] simultaneously always exists and can be 
written in the following form: 

00 00 
p f = M k U k ® M' k U' k Pin J] U\M\ ® u'Jm'J. (17) 

k=l k=l 

At the same time, the final state pf satisfies the following 
condition: 

trip/ = trap/ = 1/2. (18) 
As an example, consider the imperfect entangled state 
p = 0. VX^I + 0.12|^")<<p"| + O.78|0+)(<n (19) 

where \ip') = R(6i) ® /|00), \tp") = R{9 2 ) ® I\U) with 

0i = 7r/5,0 2 = -3tt/10 and R(9) = ( C ° S f Sm0 , Y 
' ' v ' V sin 6 — cos J 

The entanglement increase through 7 steps of iteration is 

shown in Fig. [TJ 

The remaining task is to show that, starting from the 

same state p; n , all final states satisfying Eq. (fT8)l have 

the same value for entanglement concurrence. 

Lemma 3. If state p/ satisfies Eq. (fl8|) . then state p^ = 

U®U'p f W ® U A also satisfies Eq. flU}. Here U, U' are 
any two unitary operators. 
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FIG. 1: (color online) The concurrence C versus number of 
iterations of single-qubit operation. 



This conclusion is obvious since a unity density opera- 
tor remains to be unity after any local unitary transfor- 
mation. 

Assume we have two different final states p u , p v ob- 
tained by using different processes from the same initial 
state, and they satisfy Eq. (fT8|) . Suppose 

00 00 
Pu = n M k U k ® M' k U' k pin [] uiM\®u'jMl 

k=l fe=l 
00 00 

pv = n ti' k v k ■ Pin . n vi^t ® v k w k 1 , 
fe=i fe=i 

where M k ,M.' k ,N k ,N k arc projective operators and 
U k ,U' k ,V k ,V' k are unitary operators. By using singular- 
value decomposition, We have 

Pv = wvw®w'v'w' >„-w t P t vv t ®W t 7'' t yv' t , (20) 

where VV, W, VV', W are unitary operators and V, V are 
projective operators defined in Eq. ([TJ. Denote p w = 
VV* ® VV'V^VV ® VV' and p«, = W® W'p«W+ ® W'+ . We 
have 

p w =V®V'p w V i ®V'^. (21) 

According to Lemma 3, we know that p w ,p w satisfy 
Eq. IU1). Thus "P and V' must be either identity or 
diag[l, i]. This indicates that V and T 3 ' are unitary there- 
fore the entanglement concurrence of p w and p w must be 
same. We now obtain the major result of this letter: 
Theorem 2. Given any inseparable two-qubit ini- 
tial state pi n , the entanglement concurrence can be im- 
proved through single-qubit operations provided that 
the reduced density matrix of any one qubit is not 
diag[l/2, 1/2]. Among all out-come states {p ou t|pout = 
e ® e'(pin)} through positive-definite local maps, the 
state p ou t = Q ® Q'pinQ^ ® Q has the largest en- 
tanglement concurrence if the density matrices of each 
qubit of the outcome state are 1/2. The corre- 
sponding local maps at each side are simply positive- 
definite matrices Q, Q' which can be constructed by 
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FIG. 2: (color online) The concurrence C versus 8\ and 8 2 . 
Here we have Q = M(ai, &i)I/(0i), Q' = M(a 2 , b 2 )U(8 2 ) with 
ai = 0.99987, a 2 = 0.01797, (7(6>) = 7?(6>)cr z and er z is the 
Pauli-z matrix. The peak point indicates the maximum con- 
currence 0.8858 with 9i = 0.9427 and 62 = 0.9428. 



Eq. (hzd, i.e., q = nr=i M k u k , c 



-n^i^£ where 
MkUk {M' k U' k ) diagonalize the state of the first (sec- 
ond) qubit into the form of 1/2 at the corresponding 
step. Specifically, tr 2 [(M k U k ® 7)p / fe _ 1 (M fe ?7 fe ® I) 1 ] = 
7/2 and tn[(J ® M' k U' k )p k (I <g> Mj^)*] = 1/2 where 
4 = Ilti ® MlUDp-^MiUi ® M/[//)t and w = 

M fc Z/ fc ® ///,._! (M fc C/ik)t ® 7. 

Remark: In the theorem, we have presented a mathemat- 
ical way to construct Q, Q' by iteration. We emphasize 
that, in applying our theorem in a real experiment, one 
can compute Q, Q' and then realize the physical process 
in only one step. 

Proposed experiment and numerical simulation. — We 
propose to test Theorem 2 with the initial state p as 
defined in Eq. (|19p . With many iterations, we have Q = 



n k M k u k 



0.0095 

Then we find that the entangle- 



0.0130 
0.0106 -0.0145 
0.8091 0.5874 

ment concurrence of the final state Q <8> Q!pQ} <2> Q' T is 
0.8858. Changing matrices Q and Q', the outcome entan- 
glement is always smaller than 0.8858. Numerical results 
are presented in Fig. (2) and Fig. (3). 

Concluding remark. — In summary, we have presented 
explicit results on probabilistically improving and max- 
imizing the quantum entanglement of a mixed state 
through single-qubit operations only. Testing schemes 



are proposed with numerical simulations. The local op- 
erator maximize the outcome entanglement concurrence 
and can be constructed numerically by iteration. It is 
interesting to construct the operators directly from the 
initial p- ln analytically. 
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FIG. 3: (color online) The concurrence C versus 6\ [in (a)] 
and (9 2 [in (b)]. Here we have Q = M(ai,bi)U(6i), Q' = 
M(a 2 ,b 2 )U(8 2 ) where U{0) = R{0)o z . We set 01 = 
0.99987, a 2 = 0.01797 [in (a)] and a 2 = 0.01797, 9\ = 0.9427 
[in (b)]. The peak points indicate the maximum concurrence 
0.8858 with 81 = 0.9427 [in (a)] and 8 2 = 0.9428 [in (b)]. 
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